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Abstract: We introduce the boundary measure at scale r of a compact subset of the 
n-dimensional Euclidean space. We show how it can be computed for point clouds and 
suggest these measures can be used for feature detection. The main contribution of this 
work is the proof a quantitative stability theorem for boundary measures using tools of 
convex analysis and geometric measure theory. As a corollary we obtain a stability result 
for Federer's curvature measures of a compact, allowing to compute them from point-cloud 
approximations of the compact. 
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Stabilite de mesures de bord 



Resume : Nous introduisons la notion de mesure de bord d'echelle r d'un sous-ensemble 
compact de l'espace euclidien de dimension n. Nous montrons comment calculer ces mesures 
pour un nuage de points et suggerons que ces mesures peuvent etre utilisees pour de la 
detection de features. La principale contribution de ce travail est la demonstration d'un 
theoreme quantitatif de stabilite des mesures de bord, utilisant des outils de Panalyse convexe 
et de la theorie geometrique de la mesure. En corollaire, wous obtenons un resultat de 
stabilite des mesures de courbure d'un compact (notion introduite par Federer), permettant 
de les calculer a partir d'approximations du compact par des nuages de points. 

Mots-cles : detection de dimension, nuages de points, mesures de courbure, fonctions 
convexes, plus proche voisin. 
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Introduction 

Motivations and previous work. The main goal of our work is to develop a framework 
for features detection: finding the boundaries, sharp edges, corners of a compact set K C W 1 
knowing only a possibly noisy point cloud sample of it. 

This problem has been an area of active research in computer science for some years. 
Many of the currently used methods for feature and dimension detection (see [DGGZ03] 
and the references therein) rely on the computation of a Voronoi diagram. The cost of 
this computation is exponential in the dimension and cannot be practically realized for 
an ambient dimension much greater than three. In low dimension, several methods have 
been invented for boundary detection (mostly to detect holes), for example [FK06J (2D, 
graph-based), [BSKJ (3D), and [RBBK06J. Sharp edges detection has also been studied in 
|GWM01j . and recently in |DHOS07| . 

The algorithms we develop have three main advantages: they are built on a strong math- 
ematical theory, are robust to noise and their cost depend only on the intrinsic dimension 
of the sampled compact set. None of the existing methods for feature detection share these 
three desirable properties at the same time. 

Boundary measures and their stability. Given a scale parameter r, we associate to 
each compact subset K of E™ a probability measure /?k>- This boundary measure of K at 
scale r as we call it, gives for every Borel set A C W 1 the probability that the projection on 
K of a random point at distance at most r of K lies in A (the projection on K, denoted by 
Pk, maps almost any point in R™ to its closest point in K). 

Intuitively, the measure 0k ,r will be more concentrated on the features of K: for instance, 
if K is a convex polyhedron in R 3 , (3k,t will charge the edges more than the faces, and the 
vertices even more (see example [I]). It should also be noticed that this measure is closely 
related to Federer's curvature measures (introduced in |Fed59j ). 

This article focuses on the stability properties of the boundary measures, showing that 
they can be approximated from a noisy sample of K . The problem of extracting geometric 
information from these boundary measures will be treated in an upcoming work. The main 
stability theorem can be stated as follow: 

Theorem IjlV.ip . // one endows the set of compact subsets of R™ with the HausdorfJ dis- 
tance, and the set of compactly supported probability measures on R™ with the Wasserstein 
distance, the map K i— > 0K,r is locally 1/2-Holder. 

In the sequel we will make this statement more precise by giving explicit constants. A 
very similar stability result for a generalization of Federer's curvature measures is deduced 
from this theorem. We deduce theorem IIV.1I from the two theorems IIII.5I and III. 31 below, 
which are also interesting in their own. 

Theorem i|III.5p . Let E be an open subset ofW L with (n— l)-rectifiable boundary, and f,g 
be two convex functions such that diam(V/(£') U Vg(E)) ^ k. Then there exists a constant 
C(n, E, k) depending only on n and E such that for \\f — gW^ small enough, 

l|V/-V 5 || L1(B) ^C(n,E,k)\\f~g\\ 1 J 2 
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Theorem l|II.3|) . If K is a compact set o/R™, for every positive r, dK r = {x ; d(x, K) = r} 
is (n - l)-rectifiable and H n - l {dK r ) ^ Af(dK, r) x w n _i(2r) 

Theorem IIII. 5 1 is used to show that the map K i— * px G (where px is the projection 
on if) is locally l/2-H61der, which is the main ingredient for the stability result. Theorem 
111.31 improves upon [OP85] , in which Oleksiv and Pesin prove the finiteness of the measure 
of the level sets of the distance function to K . It is used here as a tool to show that K r AK' r 
is small when K and K' are close (AAB being the symmetric difference between A and B, 
and K r being the set of points at distance at most r from K) . 

Outline. In the first section we give some examples of boundary measures and show how 
they can be computed efficiently for point clouds. The second and third sections contain the 
proofs of theorems III. 31 and IIII. 51 respectively. In the fourth section we deduce from these 
theorems the stability results for boundary and curvature measures. 

I Definition of boundary measures 

Some examples of boundary measures 

Notations. If if is a compact subset of R™, the distance to K is denned as A.k{x) = 
minyg/f \x — y|| . The r-tubular neighborhood or r-offset around a subset F C R n is the 
set of points at distance at most r from F, and is denoted by F r . 

For x £ R™, the set of points y £ K that realizes this minimum is denoted by proj^ (x). 
One can show that #projje(a;) = 1 iff d^ is differentiable at x. Since d^ is 1-Lipschitz, a 
theorem of Rademacher ensures that both conditions are true for almost every point x € W 1 . 

This allows us to define a function pk € L 1 1 oc (IR n ), called the projection on K, which maps 
(almost) every point x £ 1" to its only closest point in K. The s-dimensional Hausdorff 
measure is denoted by H s ; in particular H n coincides with the usual Lebesgue measure on 
M™. 

Definition 1,1. The r-scale boundary measure Pk., r of a compact K of R™ associates to 
any Borel set A C R™ the probability that the projection of a random point at distance less 
than r of K lies in A. 

If we denote by fix,r the pushforward of the uniform measure on K r by the projection on 
K, ie. for all Borel set ACl", mA A ) = H n (p^(A) flF), then (3 K , r = H n {K r )- Vif,r- 

Examples. 1. If C — {xi; 1 ^ i ^ N} is a «point cloud», that is a finite set of points of 
R", then (3c,r is a sum of weighted Dirac measures. Indeed, if Vorc(;ri) denotes the 
Voronoi cell of Xi, that is the set of points closer to Xi than to any other point of C, 
we have 

n 

Mc.r = J2 Hn ( VoI c(*i) H C r )5 Xi 
i=l 
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2. Let S be a unit-length segment in the plane with endpoints a and b. The set S r is the 
union of a rectangle of dimension 1 x 2r whose points projects on the segment and two 
half-disks of radius r whose points are projected on a and b. It follows that 

HS,r = 2r H 1 \ s + \r 2 8 a + ^r 2 6 b 

3. Let P be a convex solid polyhedron of ]R 3 , {e 3 } be its edges and {vk} be its vertices. 
We denote by a(e.j) the angle between the normals of the two faces containing e*, and 
by K(vk) the solid angle formed by the normal cone at Vk- Then one can see that 

HP, r = H 3 | p + ri? 2 | ap + ^r 2 a( ej ) x U 1 | e . + £ r 3 K(v k )6 Vh 

j k 

4. More generally, if if is a compact with positive reach, in the sense that there exists 
a positive r such that the projection on K is unique for any point in K r , there exist 
Borel measures ($K,i)o^i<n on ^™ such that 

n 
i=0 

where u>i is the volume of the unit sphere in These measures <&K,i are called the 

curvature measures of the compact set K and have been introduced under this form 
by Federer in [Fed5^ , generalizing existing notions in the case of convex subsets and 
compact smooth submanifolds of W 1 (Minkowski's Quermassintegral and Weyl's tube 
formula, cf. |Wey39| ). 

The second and third examples show exactly the kind of behaviour we want to exhibit 
(and so does figure II. ip : the measure 0K,r can be written as a sum of weighted Hausdorff 
measures of various dimension, concentrated on the features of K: its boundary, its edges 
and its corners. This remark together with the stability theorem for boundary measures 
shows that they are a suitable tool to be used in robust feature extraction algorithms. In 
the next paragraph we show how to compute them efficiently for point clouds. 

The boundary measure of a point cloud 

A fast method for computing the boundary measures of point clouds is of crucial importance 
for practical applications. Indeed, most real-world data, either 3D (laser scans) or higher 
dimensional is given in the form of an unstructured point cloud. Since computing the Voronoi 
diagram of a point cloud has an exponential cost in the ambient dimension, we will be using 
a probabilistic Monte-Carlo method to get an approximation of the boundary measures. In 
a very general way, if /z is an absolutely continuous measure on M n , one can compute p#cM 
as shown below. The three main steps of this algorithm (I, II, and III) are described with 
more detail in the following paragraphs. 
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Input: a point cloud C = {xi}, a measure /i 

Output: an approximation of pc#M m the form ^ k(i)8xi 

I. Choose N big enough to get a good approximation with high confidence 
while n ^ N do 

II. Choose a random point X n with probability distribution fi 

III. Finds its closest point Xi in the cloud C, add 1 to n(xi) 
end while 

return [n(xi)/N] i . 

Step I. The measure /ijv — l/^Ei<iV ^-X. where (X;) is a sequence of independent 
random variables whose law are /i is called an empirical measure. The question of whether 
(and at what speed) [In converge to \i as N grows to infinity is well-known to probabilists 
and statisticians. The results of this section are not original and can probably be improved, 
they are presented here only to give proof- of- concept bounds for N. 

Theorem 1.2 (Hoeffding's inequality). If (Yi) is a sequence of independent [0, l]-valued 
random variables whose common law v has a mean met, and Y n — 0-/N) Yli<N ^ then 

¥(\Y N - m\ > e) ^ 2exp(-27Ve 2 ) 

In particular, let's consider a family (Xi) of independent random variables distributed 
according to the law pc#M- Then, for any 1-Lipschitz function / : R™ — * K with H/H^ < 1, 
one can apply Hoeffding's inequality to the family of random variables Yi — f(Xi) : 

P 

This kind of estimate also follows from Talagrand Ti(A)-inequalities, in which case the 
factor 2 in the exponential is replaced by 2A. Bolley, Guillin and Villani use this fact to get 
quantitative concentration inequalities for empirical measures with non-compact support in 
|BGV07| . 

We now let BL 1 (C) be set of Lipschitz functions / on C whose Lipschitz constant Lip / 
is at most 1 and WfW^ ^ 1. We let A/'(BL 1 (C), H-H^ ,s) be the minimum number of balls 
of radius at most r (with respect to the H-H^ norm) needed to cover BL 1 (C). Proposition 
11.31 gives a bound for this number. It follows from the definition of the bounded-Lipschitz 
distance (see II. 4|) and from the union bound that 

P[d bL (pc#MJV,Pc#M) >e\ ^2AA(BL 1 (C7),||.|| 0O ,e/4)exp(-7V £ 2 /2) 
Proposition 1.3. For any compact metric space K , 

^(BLHmui^x f-J 



N 

*E /(*>-/ 



fdfi 



< 2exp(-2iVe 2 ) 
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Proof. Let X = {xi} be an e/4-dense family of points of K with #X = J\f(K,e/i). It is 
easily seen that for every 1-Lipschitz functions /, g on K, \\f — gW^ ^ || (/ — fiOlxlloo + 
Then, one concludes using that N(Bl}(X), IUU ,e/2) < (A/e)* x . □ 

In fine one gets the following estimate on the bounded-Lipschitz distance between the 
empirical and the real measure: 

P[d bL (pc#Mjv,Pc#m) >e] ^2exp(ln(16/e)AA(C7,e/16)-7V£ 2 /2) 

Since C is a point cloud, the coarsest possible bound on 7V(C, e/16), namely #C, shows 
that computing an e-approximation of the measure p#^i with high confidence (eg. 99%) can 
be done with N = 0(#C ln(l/e)/e 2 ). 

Step II. To simulate the uniform measure on K r one cannot simply shoot points in a 
bounding box of K r , keeping those that are actually in K r since this has an exponential cost 
in the ambient dimension. Luckily there is a simple algorithm to generate points according 
to this law which relies on picking a random point x% in the cloud C and then a point X in 
B(xi, r) — taking into account the overlap of the balls B(x,r) where x E C: 

Input: a point cloud C — {xi}, a scalar r 
Output: a random point in C r whose law is H n \ Kr 
repeat 

Pick a random point Xj in the point cloud C 
Pick a random point X in the ball B(xi,r) 

Count the number k of points Xj € C at distance at most r from X 

Pick a random integer d between 1 and k 
until d = 1 
return X. 

Step III. The trivial algorithm for computing the projection of a point on a point cloud 
takes exactly n steps. Since generally N will an order of magnitude greater than n we 
might improve the overall 0(n 2 ) cost by maintaining a data structure which allows fast 
nearest-neighbour queries. This problem is notoriously difficult and until recently most of 
the efficient algorithms in high dimension were only able to compute approximate nearest 
neighbours. This amounts to replacing pc by a map p e with the property that for all x, 
||p e (x) — pc(x)|| ^ (1 + e)dc(x). Unfortunately, the techniques we develop in this paper 
do not seem to apply directly to get quantitative closeness estimates for the measures p e #/U 
and Pit#^- 

It should be noted that for low entropy point clouds, nearest neighbor queries can be done 
more efficiently. For instance, a recent article by Beygelzimer, Kakade and Langford (cf. 
|BKL06| ) introduces a structure called cover trees which allows an exact nearest neighbour 
query with complexity 0(c 12 logn) where c is related to the intrinsic dimension of the point 
cloud, with an initialisation cost of 0(c e nlogn). 
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Figure LI: Boundary measure for a sampled mechanical part. 



Wasserstein distance and stability 

Since our goal is to give a quantitative stability result for boundary measures, we need 
to put a metric on the space of probability measures on K n . The Wasserstein distance, 
related to the Monge-Kantorovich optimal transportation problem seemed intuitively (and 
later happened to really be) appropriate for our purposes. A good reference on this topic is 
Cedric Villani's book |Vil03j . 

Definition 1.4. The set of measures (resp. probability measures) on R" is denoted by 
M(R n ) (resp. M 1 (R n )). We endow M(S. n ) with the bounded Lipschitz distance, ie. 

yfi,is £ M(W l ),d hL ((i.is) = sup 

llvl| Lip ^i 

where the supremum is taken over all Lipschitz functions ip with |M| Lip = Lip<y3 + Hyll^ ^ 1 
(Lipyp being the smallest constant k such that ip is /c-Lipschitz). 

We put two distances on 7W 1 (M") (which are in fact identic, see below). The Fortet- 
Mourier distance, which is almost the same as the bounded Lipschitz one: 

V/i, v G A4 1 (M. n ), c^fmCMi v ) ~ sup 

Lip ip ^ 1 

And the Wasserstein distance: 

Wi(jj.,u) = inf {E(d(X,Y)) ; law(A) = fi, law(Y) = v] 

where the infimum is taken over all random variables X and Y whose laws are fi and v 
respectively. 



ipdfj, — / (fidis 



ipdfi — / (pdv 
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Notations. If fi and v G A4(M. n ) are absolutely continuous with respect to 7i n , ie. d/i = 
ipdH n and dv = ipdTL n we denote by fiCw the measure defined by d(^n^) = min((^, ip)dTL n , 
and fiAv = /x + v — 2/j, fl v. 

Proposition 1.5. If ^ & A^(IR™) is absolutely continuous with respect to the Lebesgue mea- 
sure, and f,g : 1™ — ► R n are two functions in then 

// /i and f are iwo absolutely continuous measures on W 1 , 

d hL (.f#(J-,g#v) ^ ||/-5||i( M rv) + mass(/zAi/) 
Proof. For any 1-Lipschitz function <p on R", 

J <pdf # IJ,- J (fdg # ii = J ipofdfl- J (fiogdfi 

Lip^ J ||/-5l|dM<||/-5l|i M 



For the second inequality, let us first remark that there exists two positive measures \i r 
and v r such that jj, = n fl v + fj, r and v = /i n v + v r . Then, 

< d bL (/#M, /#M n i/) + d hL (f # [i n i/, 3 # M H z/) + d hL (g#n, g # ^ n i>) 
Now let us bound one of the extreme terms of the sum, 



v^s.tiML ^ i, 



J (fdf # n- J ipdf # iiC\v 



ip O fd/Ir 



mass(/x r ) 



One concludes using that \x r + f r = /xAz/. 

Corollary 1.6. //if and K' are two compact subsets ofW L , 

4L(MKr,^,r) < \\p K - Vk II i { K^K'r) + H n (K r AK' T ) 



□ 



Hence to get a quantitative continuity estimate for the map K i— > /iif ir one needs to 
show that if if and if' are Hausdorff-close, K r AK' r is small, and to evaluate the continuity 
modulus of K \— » p^ e 1 (if r n if /r ). This is the purpose of the two following paragraphs. 



RR n° 6219 



10 



Chazal, Cohen-Steiner & Merigot 



II K r AK' r is small when K and K' are close 

It is not hard to see that if dn{K, K') is smaller than e, then K r AK' r is contained in 
(K r+£ \ K r ~ £ ). The volume of this thick tube around K can then be expressed as an 
integral of the area of the hypersurfaces dK* . 

The next proposition gives a bound for the measure of the r-level set dK r of a compact 
set K C R™ depending only on its covering number Af(K, r) {ie. the minimal number of 
closed balls of radius r needed to cover K). In what follows, K r is the set of points of R" 
at distance less than r of K, and dK r is the boundary of this set, ie. the r-level set of &k- 
In this paragraph, we prove the following theorem : 

Theorem. If K is a compact set o/R™, for every positive r, dK r is 7i" _1 -rectifiable and 
H n -\dK r ) ^JV(dK,r) x w n _i(2r) 

This proposition improves over a result of finiteness of the level sets of the distance 
function to a compact set, proved by by Oleksiv and Pesin in [OP85J . We begin by proving 
it in the special case of "r-flowers". A r-flower F is the the boundary of the r-tube of a 
compact set contained in a ball B(x,r), ie. F = dK r where K C B(x,r). The difference 
with the general case is that if K C B(x, r), then K r is a star-shaped set with respect to x. 
Thus we can define a ray-shooting application sk ■ S n ^ 1 — > dK r which maps any v e S n ^ 1 
to the intersection of the ray emanating from x with direction v with dK r . 




Figure II. 2: Ray-shooting from the center of a flower. 



Lemma II. 1. Let K = {e} C B(x,r) and define s e as above. Then s e is 2r-Lipschitz (with 
respect to the sphere's inner metric) and its Jacobian is at most (2r)™ _1 . 

Proof. Solving the equation + tv — e\\ — r with t gives 

s e (v) = x + (v\x — e) 2 + r 2 — \\x — e|| 2 — (v\x — e) \ v 

Denote by H v the orthogonal of the 2-plane P spanned by v and s e (v) — e. For each 
vector w chosen in H v , a simple calculation gives: 

s e (v + tw) = s e (v) + tw \\s e (v) — x\\ + o(t 2 ) 
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Hence the derivative of s e along H v is simply the multiplication by ||s e (w) — x\\ ^ 2r. 

Now, we now consider the case of the 2-plane P. We denote by 9 the angle between 
s e (v) — x and s e (v) — e and by w a vector tangent to v in the intersection of the sphere with 
P. Then 

\\(ds e ) v (w)\\ _ \\s e (v) - x\\ 
\\w\\ \co S (0)\ 

Now let us remark that 

||« e («) - e|| \\s e (v) - x\\ \cos(0)\ = \(s e (v) - e\s e (v) - x)\ 

= \{\\x-s e {v)\\ 2 + \\se{v)-e\\ 2 -\\x-e\\ 2 ) 
>\\\x-*Mt 

Finally we have proved that ||(ds e )„|| ^ 2r. The result follows by integration. □ 

We denote by uj n (r) the n-Hausdorff measure of the n-sphere of radius r. 

Corollary II. 2. A r -flower in R™ is a 7i™ _1 -rectifiable set and its measure is at most 
w„_i(2r). 

Proof. Let K C B(x,r) be the compact set generating the flower dK r . As above, for any 
vector v G S"^ 1 , we denote by s the intersection of the ray {x + tv ; t > 0} with dK r . Since 
if r is a star-shaped set around x, s is a bijection from S" -1 to dK r . 

Now let (yk) be a dense sequence in K, and denote by Sk the projection from <S™ -1 to 
the flower 9(Ui^fc{j/i}) r defined as above. Then (sk) converges simply top on Indeed, 
if we fix v € and e > 0, the segment joining x and s(v) truncated at a distance e of 

s(v) is a compact set contained in mtK r . It is covered by the union UjS(t/j,r), so that for 
N big enough it is also covered by Uk^NB(yk,r). For those N, \\sk(x) — s(x)\\ < e. 

Finally, dK r is the image of the sphere by p, which is 2r-Lipschitz as a simple limit of 
2r-Lipschitz functions. □ 

We now deduce a general bound on the measure of the tube boundary dK r around a 
general compact set K by covering it with a family of flowers: 

Theorem II. 3. If K is a compact set ofW 1 , for every positive r, dK r is a Ti. 71 ^ 1 -rectifiable 
subset ofW 1 and moreover, 

H n -\dK r ) < N{dK ir ) x w„_i(2r) 

Proof. It is easy to see that dK r C d{dK r ). Thus, if we let (xi) be an optimal covering 
of dK by open balls of radius r, and denote by Ki the (compact) intersection of dK with 
B(xi,r), the boundary dK r is contained in the union V3idK\. Hence its Hausdorff measure 
does not exceed the sum J2 i TL n ^ 1 {dKl). One concludes by applying the preceding lemma. 

□ 
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Remark II. 4. 1. The bound in the theorem is tight, as one can check taking K = B(0, r). 

2. Let us notice that for some constant C(n), J\f(B(0, l),r) ^ 1 + C(n)r~ n . From this 
and the above bound it follows that 

H n ~ 1 (dK r ) ^ (1 + C(n) x (diam(X)/r)"K_ 1 (2r) 

^C»x(l + diam( ^ )n ) 
r 

for some universal constant C'(n) depending only on the ambient dimension n. This 
last inequality was the one proved in [OP85J. 

To conclude we use a weak formulation of the co-area formula, a standard result of 
geometric measure theory (|DG54|, [Fcd59|), which reads 

J \V x f\ dH n (x) = ^H"- 1 (r 1 ( 2 /))dH 1 (y) 

whenever / : R™ — > R is a Lipschitz map. From this formula and the previous estimation 
follows that 

Corollary II. 5. For any compact sets K,K' C R™, with dii(K,K') ^ e, 
H n (K r AK' r ) sC I + H'^idK^dt 

J r — e 

< 2N{K, r - e)cj„_i(2r + 2s) x e 

III The map K i— > p% is locally l/2-H61der 

We now study the continuity modulus of the map K t— > p K e 1 (i?) , where E is a suitable 
open set. We remind the reader of two well-known facts of convex analysis (see for instance 
pa83]): 

1. If / : SI C R™ — > R is a locally convex function, its subdifferential at a point x, 
denoted by d x f is the set of vectors v of R" such that for all h e R" small enough, 
f(x + h) ^ f(x) + (h\v). Then / admits a derivative at x iff d x f — {v} is a singleton, 
in which case V x / = v. 

2. A locally convex function has a derivative almost everywhere. 

Lemma III. 1 . The function vk ■ R™ — > R, x ^ \\x\\ 2 — dnix) 2 is convex with gradient 
Vvk = %Pk almost everywhere. 

Proof. By definition, v K (x) = sup yeK \\x\\ 2 - ||a;-?;|| 2 = sup^^ v K , y {x) with v K , y (x) 
' — \\y\\ 2 . Hence vk is convex as a supremum of affine functions. Because vx. PK ( x ) 
and Vk take the same value at x, d x VK. PK ( x ) = {%Pk (x)} C 8vk ■ Since vk is differentiable 
almost everywhere, equality must be true almost everywhere which concludes the proof. □ 
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This lemma shows that | \pk ~ Pk' || 1 = 1/2 || Vvr- — Vwk' | j i , E s . Our estimation of the 
continuity modulus of the map K will follow from a general theorem which asserts 

that if <p and ^ are two uniformly close convex functions with bounded gradients then V<p 
and V-0 are 1 -close. The next proposition below is the 1-dimensional version of this result, 
from which we then deduce the general theorem. 

Proposition III. 2. If I is an interval, and ip : I — > K and tp : / — > M are iiuo convex 
functions such that di&m(ip' (I) U ip'(I)) ^ fc, iften letting 5 = \\(p — ip\\oa(i\, 

J \tff - f//\ < 67r(length(7) + fc + (J 1 / 2 )^ 2 

Lemma III. 3. Let / : I — » R be a nondecreasing function with diam </?(/) ^ fc. Then, if F is 
the completed graph of /, ie. the set of points (x,y) € JxR such that lim.^- ip ^ y ^ lim x + 93, 
then H n (F r ) < 37r(length(7) +Hr)xr. 

Proof. Let 7 : [0, 1] — > i* 1 be a continuous parametrization of F, increasing with respect to the 
lexicographic order on M 2 . Then, for any increasing sequence (ti) G [0, 1] and (a;,, j/i) = 7(^), 

X! IIt(*»+i) - 7(*<)ll < _a: * + y-'+i _ < length(J) + fc 

Hence length (F) ^ length(J) + fc. Thus we can choose a 1-Lipschitz parametrization of 
F, 7 : [0, length(J) + k] — > F. Then for any positive r, the set X = xr);0^K ^V} 
with N the upper integer part of (length(J) + k)/r, is such that any point of F is at distance 
at most r of X. Hence F r is contained in X 2r , implying that H n {F r ) 7V7r(3r/2) 2 
37r(length(7) + fc + r)r. □ 

Proof of proposition UIT7B . Let I — [a, 6] and J = [c, c + fc] be such that <//(/) U V>'(-0 £ J. 
Without loss of generality we will suppose that t/j'(a) — ip'(a) — c and ip'(b) = <p'(b) = c + k. 
With this assumption, the completed graphs <£> and \P of ip' and ip' defined as above are two 
rectifiable curves joining (a, c) and (6, c+ k). We let F be the set of points (x, y) € K 2 lying 
between those graphs; the quantity we want to bound is L \ip' — tp'\ = 7i 2 {V). 

Let 8 = \\ip — ip\\°v(jy For any point p — [x,y) in V, and any 8' > 8, the closed disk 

D = B(p, y/28 1 /-k) of volume 28' centered at p cannot be contained in V. Indeed if it were, 
then the difference k — ip — ip would increase too much around p: since k' has a constant 
sign on this segment, 

\k(x + 28' jn) - k(x - 28' jix)\ = / \n'\ ^ H 2 (D) = 28' > 28 

Jx-15' fit 

This contradicts H^H^ = 8. Hence, D must intersects dV implying that V must be contained 
m for any 8' > 8. Since dV = $U$, the previous lemma gives 

H 2 (V) sj H 2 (^V 2 ^^ + n 2 (^V 2 ^) 67r(length(7) + k + y/WJj^y/WJw 

Letting 8' converge to 8 concludes the proof. □ 
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A generalization of this proposition in arbitrary dimension will follow from an argument 
coming from integral geometry, ie. we will integrate the inequality of proposition IIII.2I over 
the set of lines of R™ to get a bound on \\\/<p — V^||i( B ). 

We let L n be the set of oriented affine lines in R™ seen as the submanifold of R 2 ™ made 
of points (u,p) £ R™ x 1™ with u <E <S n_1 and x in the hyperplane {u} ± , and endowed with 
the induced Riemannian metric. The corresponding measure dC n is invariant under rigid 
motions. We let 2?" be the set of oriented lines with a fixed direction u. 

The usual Crofton formula (cf. [Mor88j for instance) states that for any H n_1 -rectifiable 
subset S of R™, with /?„ the volume of the unit n-ball, 

H n -\S) = — ^— [ #{t^S)d£ (III.l) 

where #X is the cardinality of X. We will also use the following Crofton-like formula: if K 
is a W n -rectiflable subset of R", 

H n (K) = / 7i 1 (£nX)d£ (III.2) 

which follows from the Fubini theorem (remember u> n -i is the volume of the (n— l)-sphere). 
Lemma III. 4. Let X : E — > R™ be a 1 -vector field on an open subset E C 1". 

/ pr|| = -JL_ / / \{X{y)\u{i))\dyM 

JE j££C n Jy££nE 

Sketch of proof. The family of vector fields of the form J^i ^iXs^i j where the fli are a finite 
number of disjoint open subsets of R™ and Xi are constant vectors, is 1 -dense in the space 
1 (R n ,R n ). Using this fact and the continuity of the two sides of the equality, it is enough 
to prove this equality for X — x \\X\\ \e where x is a constant unit vector and E a bounded 
open set of R" . 

In that case, one has 

f f \(X(y)\u)\dyd£=\\X\\\(x\u)\ [ length^ n l)dt 

= \\X\\ 1{E) \(x\u)\ 

By a Fubini-like theorem one has 

/ f \(X(y)\u(e))\dyde= [ [ [ \{X(y)\u(£))\dydedu 

= \\X\U( E )[ \(x\u)\du 
'Jugs™- 1 
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The last integral does, in fact, not depend on x and its value can be easily computed: 



f |(x|u)|du = 2cj„_2 f til-t 2 )^- 1 ^ 



2 

— U n -2 



□ 



Theorem III. 5. Let E be an open subset ofW 1 with (n — l)-rectifiable boundary, and f,g 
be two locally convex functions on E such that diam(V/(i?) U Vg(E)) ^ k. Then, letting 
& = II/- fflU(ij) 

II V/ - V.g|| Ll(£) < Ci(n)(« n (£) + (k + S^TT-^dE)^ 1 ' 2 

with Ci(n) ^ Gnn as soon as n > 5 (in fact, C\(n) = 0{\fn) ). 

Proof of the theorem. The 1-dimensional case follows from proposition IIII.2I in that case, 
E is a countable union of intervals on which / and g satisfy exactly the hypothesis of the 
proposition. Summing the inequalities gives the result with Ci(l) = 67r. 

The general case will follow from this one with the use of integral geometry. If we set 
X = V/ - Vg, fe = f\ enE and g e = g\ enE . Lemma HE! gives, letting D(n) = n/(2w„_ 2 ), 

[ \\Vf-Vg\\ = D(n) [ [ |(V/ - Vg\u(£))\ dyd£ 

JE J£GC n JyGinE 



= D(n) / \ti-g' e \dyd£ 

Ji£C n JyGinE 

The functions ft and gi satisfy the hypothesis of the one-dimensional case, so that for 
each choice of t, and with 8 = \\f — g\\ao/ E \, 

|/; - g'f\ dy < &xD(n)(H\E n£) + (k + 5 l ' 2 )U\dE n e))5^ 2 

yeinE 

It follows by integration on C n that 

J ||V/ - V5|| < 67rD(n) ^ H^E n £)dC n + (k + S 1/2 ) J H°{dE n £)d£ n ^j 5 1/2 

The formula IIII.ll and IIII.2I show that the first integral is equal (up to a constant) to the 
volume of E and the second to the (n — l)-measure of dE. This proves the theorem with 
C\(n) = 6irD(n)(u n -i + 2/3 n _i). To get the bound on C\{n) one uses the formula w n — 1 = 
n(i n and /3 n +i ^ n as soon as n > 5. □ 

Multiplying / and g by the same positive factor t and optimizing the result in t yields a 
better, homogeneous, bound : 
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Corollary III. 6. Under the same hypothesis as in theorem I III. 51 one gets the following 
bound, with 5 = \\f — g\\aotE) ■' 

IIV/ - Vg|| L1(B) < 2C 1 (n)[(H n (E)n n - 1 (dE) diam(V/(£) U Vg^))) 1 / 2 

+ H n - 1 {dE)5 1/2 ]5 1/2 

Remark III. 7. To get an homogeneous bound as in this corollary, one could also optimize 
the one-dimensional bound of proposition IIII. 21 before integrating on the set of afHne lines of 
K" as in the proof of theorem 1111.51 The bound obtained this way is always strictly better 
than the ones of both theorem IIII. 51 and corollary IIII.6| but involves an integral term 

f ^JH.°^ndE)n 1 {lr\E)At 

whose intuitive meaning is not quite clear. 

Applying theorem IIII. 51 to the functions vk and vk 1 introduced at the begining of this 
part and using lemma [Til.!} one easily gets : 

Corollary III. 8. If E is an open set ofM. n with rectifiable boundary, K and K' two compact 
subsets ofW 1 then, with Rk = ||oo:(gj and e = dn{K, K'), 

Wpk - PK'\U(E) < Ci(n) [H n (E) + (diam(tf) + e + (2R K + s^e^H^idE)} 

x (2R K +e) 1 / 2 e 1 / 2 

In particular, if d# (K, K') is smaller than mm^Rx, diam(i^), dia.m(K) 2 / Rk ), there is an- 
other constant C2(n) depending only on n such that 

I \pk -PK'h {E) < C 2 (n) [H n (E) + dmm(K)H n - 1 (dE)} ^R K A H {K,K') 

Remarks III. 9. 1. This theorem gives in particular a quantitative version of the conti- 
nuity theorem 4.13 of |Fed59| : if (K n ) is a sequence of compact subsets of R n with 
reach(i4T„) ^ r > 0, converging to a compact set K, then reach(_ftT) > r and pk„ 
converges to pk uniformly on each compact set contained in {x € K n ; &k{x) < r}. 
However we have to stress that the result we have proved is more general since it 
does not make any assumption on the regularity of K n — at the expense of uniform 
convergence. 

2. The second term of the bound involving H. n ~ 1 (dE) is necessary. Indeed, let us suppose 
that a bound \\pk — Pk'\\i(e) ^ C(K)H. n (E)y/e were true around K for any open set 
E. Now let K be the union of two parallel hyperplane at distance R intersected with 
a big sphere centered at a point x of their medial hyperplane M . Let E e be a ball of 
radius e tangent to M at x and K e be the translation by e of K along the common 
normal of the hyperplanes such that the medial hyperplane of K e touches the ball E E 
on the opposite of x. Then, for e small enough, \\pk — Pk'\\i(e c ) — R x which 
clearly exceeds the assumed bound for a small enough e. 
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3. According to this theorem, the map K i-> p K e 1 (E) is locally l/2-H61der. The 
following example shows that this result cannot be improved even around a very simple 
compact set. 




R 



Figure III. 3: A sequence of «knife blades» converging to a segment. 

Let S and S' be two opposite sides of a rectangle E, ie. two segments of length L and 
at distance R. We now define a Hausdorff approximation of S: for any positive integer 
N, divide S in N small segments s* of common length £, and let Cj be the unique 
circle with center in S' which contains the two endpoints of s l . We now let Sn be the 
union of the circle arcs of Cj comprised between the two endpoints of s l . 

Then it is not very hard to see that if R E = R + e is the common radius of all the C % , 
Rl = R 2 + (£/2) 2 , ie. d H (S,S N ) = y/R 2 + (£/2) 2 - R ^ R£ 2 /8. Then the ^distance 
between the projections on S and Sn is at least Q(t) (because almost half of the points 
in E projects on the corners of Sn, see the shaded area in fig. 1111.3(1 . Hence, 

Wps -ps* II = m = n(d*(s, ^) 1/2 ) 

Replacing 1 (E) with where /i has bounded variation 

As we have seen before, a corollary of the previous result is that if [i = Ti. n \ E , the map 
K I— > pk#H is locally l/2-H61der. This result can be generalized when // = mW" where 
u G Lj L oc (M") has bounded variation. We recall some facts about the theory of functions with 
bounded variation, taken from |AFP00| . If fi C M" is an open set and u S L 1 1 oc (fi), the 
variation of u in fl is 

v(«,n) = sup|jf udiv^^GC^n),!!^ < i| 

A function u 6 L 1 1 oc (ri) has bounded variation if Q) < +oo. The set of functions of 
bounded variation on fl is denoted by BV(57). We also mention that if u is Lipschitz on fl, 
then V[u,0) = \\ Vu||i^. Finally, we let V(u) be the total variation of u in R". 

Theorem 111.10. Lei /i e M(R n ) 6e a measure with density u € BV(R") with respect to 
the Lebesgue measure, and K be a compact subset o/R™. We suppose that supp(u) C K R . 
Then, if dn(K, K') is small enough, 

d hL (PK#V,PK>#n) < C 2 (n) (\\u\\ 1(KR) +di&m{K)Y{u))VRy.A H {K,K')^ 2 
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Proof. We begin with the additional assumption that u has class C°° . The function u 
can be written as an integral over t € M of the characteristic functions of its superlevel sets 
E t = {u > t}, ie. u(x) = J °° XE t (x)dt. Fubini's theorem then ensures that for any Lipschitz 
function / defined on K™ with ||/|| Lip *S 1, 

PK'#/x(/) = / f op K ,(x)u(x)dx 



f ° Pk ' (z)X{ u >t} (x)dxdi 

By Sard's theorem, for almost any t, dE t — u~ l {t) is a (n — l)-rectifiable subset 
of K". Thus, for those i the previous corollary implies, for e — dH(K,K') ^ eo ~ 
min(-R, dia,m(K),dia,m(K) 2 /Rk), 

J 1/ °Pk(x) - f°1?K'(x)\ dx < ||pk - PK'\\i {Et) 

sC C 2 {n)[n n {E t )+dia,ia{K)K n - 1 {dE t )]VRl 
Putting this inequality into the last equality gives 

\ P K#Kf) - PK'Mf)\ < C 2 {n) ^H n (E t ) + di a m(K)H n ~\dE t )d?J Vite 
Using Fubini's theorem again and the coarea formula one finally gets that 

\pK#v{f) ~PK'#fJ-{f)\ < c 2(n) + diam(if) V(u)) ^/Re. 

This proves the theorem in the case of Lipschitz functions. To conclude the proof in the 
general case, one has to approximate the bounded variation function u by a sequence of C°° 
functions (u n ) such that both \\u — u n \\i, KR \ and \V(u) — V(u„)| converge to zero, which is 
possible by theorem 3.9 in [AFPOOj . □ 



Remark III. 1 1 . Taking u — xe where E is a suitable open set shows that theorem 1111,81 can 
also be recovered from IIII.101 



IV Stability of boundary and curvature measures 

We combine the results of corollaries 11.61 IH.5I and IIII.8I to get 

Theorem IV. 1. If K and K' are two compact sets with e = dn(K,K') smaller than 
min(diam K, r, r 2 / diam K), then 

d h h{^K.r,^K',r) < C 3 {ri)N{K, r - e)r n [r + diam(iT)] 

In particular, if for a given bounded Lipschitz function f on R n , one defines ifK,f(r) — 
l*K,r{f)> the map K i-> ip K j e C°([r min , r max ]) with < r min < r max is locally 1/2-Holder. 
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In what follows we suppose that (rj) is a sequence of n distinct numbers < rg < 



< r v 



For any compact set K and / <G 
system 



"), we let 



Vi s.t < i < n, ^ W „_ j ^ ) .(/)rr J - Mat ,T-j (/) 

Since the system is linear in (fJ,K,n(f)) an d these values depends continuously on /, the map 
/ i > i s a l so linear and continuous, «e. $^ is a signed measure on R™. It is also to 

be noticed that if K has positive reach with reach(if ) > r n , the coincide with the usual 
curvature measures of K. In that case, the following result gives a way to approximate the 
(usual) curvature measures of K from a Hausdorff-approximation of it even if its reach is 
arbitrary small. 

Corollary IV. 2. There exist a constant C depending on K and (r) such that for any 
compact subset K' of W 1 close enough to K , 



be the solutions of the linear 



Vi, 



4l ^cd H (K,K>y /2 
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